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Abstract 


Micromechanic  theory  of  fatigue  crack  initiation  by  Lin  and  his  associates  is  reviewed. 
Intrusions  and  extrusions  have  been  observed  in  fatigue  specimens.  An  initial  stress  field 
favorable  for  the  growth  of  extrusion  or  intrusion  can  be  produced  by  arrays  of  dislocation 
dipoles  as  shown  by  Lin  (1969).  These  dipoles  cause  an  initial  tensile  strain  giving  an  elon¬ 
gation  which  is  called  the  static  extrusion.  It  has  been  suggested  that  after  extrusion  reaches 
the  height  of  static  extrusion,  extrusion  growth  terminates.  Along  the  direction  of  an  extru¬ 
sion,  tensile  strain  and  tensile  stress  occurs.  This  tensile  stress  causes  the  resolved  shear 
stress  to  reach  critical  and  slide  in  a  secondary  slip  system.  The  slip  causes  plastic  tensile 
strain  which  increases  significantly  the  extrusion  growth.  Similarly  the  same  mechanism 
exists  for  the  growth  of  intrusions.  This  explains  why  extrusion  grows  much  more  than  the 
static  extrusion. 
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Introduction 


McCommon  &  Rosenberg  [1]  and  MacCrone  et  al.  [2]  show  that  metal  are  subject  to 
fatigue  at  temperature  as  low  as  1.7%.  This  seems  to  indicate  that  although  surface  corrosion, 
gas  adsorption,  gas  diffusion  into  the  metal  or  vacancy  diffusion  to  form  voids  can  have 
important  effect  on  but  are  not  necessary  to  the  fatigue  crack  initiation.  Local  plastic  defor¬ 
mation  resulting  from  the  displacement  of  dislocation  generally  proceeds  this  initiation. 

Forsyth  and  Stubbington  1954  [3]  reported  the  detection  of  extrusion  in  slip  bands  dur¬ 
ing  fatigue  of  some  aluminum  alloys.  Thompson,  Wadsworth  and  Louat  [4]  and  Hull  [5] 
detected  the  extrusion  process  in  both  copper  and  aluminum.  This  initiation  of  extrusion  pro¬ 
cess  was  also  observed  by  Meke  and  Blochwitz  [6]  and  Mughrabi  [7]  in  their  studies  of  per¬ 
sistent  slip  bands. 

Following  the  clue,  which  the  observation  on  extrusions  and  intrusions  in  slip  bands 
have  provided,  a  number  of  theories  of  fatigue  crack  initiation  have  been  proposed  by 
different  distinguished  investigators.  For  example,  Mott  [8]  proposed  that  a  screw  dislocation 
repeats  its  path  through  cross  slip.  He  considered  a  column  of  metal  containing  a  single 
screw  dislocation  intersecting  a  free  surface.  When  this  dislocation  travels  a  complete  circuit, 
the  volume  contained  in  the  circuit  is  translated  parallel  to  the  dislocation.  This  causes  the 
metal  to  extrude.  This  mechanism  does  not  explain  why  the  dislocation  under  cyclic  stress¬ 
ing,  does  not  oscillate  back  and  forth  along  the  same  path  rather  than  traversing  a  closed  cir¬ 
cuit.  Clearly  some  form  of  gating  mechanism  is  required  to  convert  the  back  and  forth  oscil¬ 
lations  of  screw  dislocations  into  unidirectional  circuits.  Cottrell  and  Hull  [9]  proposed  that 
Frank  Read  sources  exist  on  two  intersecting  slip  planes  and  a  complete  cycle  of  forward  and 
reversed  loading  results  in  an  extrusion  and  intrusion.  Such  a  model  would  predict  the  extru¬ 
sion  and  intrusion  to  form  in  neighboring  slip  bands  and  to  be  inclined  to  each  other,  but  they 
have  been  found  to  occur  in  the  same  slip  band  and  to  be  parallel  to  each  other.  Thompson 
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[10]  proposed  an  edge-screw  interaction  model  for  the  initiation  of  extrusion.  This  model 
assumes  a  change  of  spacing  of  a  pair  of  parallel  screw  dislocations  after  being  cut  an  edge 
dislocation.  This  cut  causes  jogs  in  both  edge  and  the  screw  dislocations  and  this  cut  will 
inhibit  the  change  of  the  spacing  of  the  pair  of  screw  dislocations  [11].  It  is  also  hard  to  see 
how  intrusions  are  formed  by  such  a  model.  McEvily  and  Machlin  [12]  proposed  a  model,  in 
which  two  screw  dislocations  terminate  in  a  surface  and  intersect  a  node  where  three  disloca¬ 
tions  meet.  Under  an  alternating  shear  stress,  the  two  screw  dislocations  are  assumed  to  shift 
around  a  circuit  causing  an  intrusion  and  an  extrusion  formed  in  the  same  slip  band.  How¬ 
ever,  this  model  as  pointed  out  by  Kennedy  does  not  explain  why  these  two  dislocations 
travel  around  a  circuit  instead  of  back  and  forth.  Wood  [13]  proposed  a  simple  model  of  a 
single  operative  slip  system.  An  undirectional  stressing  causes  layers  of  metal  to  slide  in  the 
same  direction;  but  forward  and  reverse  stressing  causes  different  amounts  of  net  slip  on 
different  planes  and  results  in  peaks  and  valleys.  However,  this  model  and  other  models  of 
random  slip  do  not  explain  why,  under  an  alternate  loading,  the  slip  continues  to  monotoni- 
cally  deepen  the  valley  and  raise  the  peaks  as  observed  in  experiments.  Drawbacks  of  other 
theories  have  also  been  discussed  by  Kennedy  [11]. 

For  a  dislocation  to  glide,  it  (1)  must  glide  along  a  certain  direction  on  a  crystal  plane 
(2)  must  subject  to  a  resolved  shear  stress  equal  or  greater  than  the  critical  shear  stress.  The 
above  mentioned  theories  show  the  possible  paths  of  dislocation  movement  to  satisfy  condi¬ 
tion  (1)  but  the  resolved  shear  stress  field  caused  by  the  dislocation  movement  that  has 
significant  effect  on  (2)  was  not  considered.  In  the  micromechanic  theory  developed  by  Lin 
and  his  associates,  [14-16],  this  important  effect  of  this  stress  field  in  causing  the  initiation  of 
extrusion  and  intrusion  is  considered. 
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Dependency  of  slip  on  resolved  shear  stress 


Single  ciystal  tests  [17]  have  shown  that  under  stress,  slip  occurs  along  certain  direc¬ 
tions  on  certain  crystal  planes.  These  directions  and  planes  are  generally  those  of  maximum 
atomic  density.  Slip  has  been  found  to  depend  on  the  resolved  shear  stress  and  not  on  the 
normal  stress  on  the  slidmg  plane.  This  dependency  of  slip  on  the  resolved  shear  stress  under 
monotonic  loadings  has  also  been  found  to  hold  under  cyclic  loadings  [18].  Hence  to  find  the 
slip  in  a  slip  band  or  fatigue  band,  we  need  to  determine  the  resolved  shear  stress  distribution 
in  the  metal.  To  determine  this  stress  field,  the  following  analogy  between  plastic  strain  and 
applied  force  is  used. 

Analogy  between  plastic  strain  and  applied  force 


The  resolved  shear  stress  to  cause  sliding  is  called  the  critical  shear  stress.  When  this 
critical  shear  stress  is  reached  in  some  region  in  a  body,  slip  occurs  and  causes  plastic  strain. 
If  the  load  is  removed,  plastic  strain  remains  and  cause  a  residual  stress  field.  To  find  this 
stress  field,  the  analogy  between  plastic  strain  and  applied  force  developed  by  Lin  [19,20]  is 
used.  This  analogy  is  briefly  reviewed  as  follows: 


Referring  to  a  set  of  rectangular  coordinates,  the  strain  component  is  composed  of  the 
elastic  part  denoted  by  single  prime  and  the  inelastic  part  denoted  by  double  prime 

eij  =  e'ij  +  e"ij 

(1) 

Thermal,  creep  and  plastic  strains  are  considered  to  be  the  inelastic  strain.  Neglecting  the 
anisotropy  of  the  elastic  constants,  the  stress  is  related  to  the  elastic  strain  as 


=  5 ,y  XQ'  +  2\le'ij, 

=  5  ft  (9  -  0')  +  2p  (<-  -  e'ij), 


(2) 
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where  X,  and  | x  are  Lame’s  constants,  0  is  the  dilation  and  9  is  the  inelastic  dilation.  The 
condition  for  equilibrium  within  a  body  of  volume  v  is 

xijj  +  Fi  =  0  in  v> 

(3) 

where  the  subscript  after  the  comma  denotes  differentiation,  the  repetition  of  the  subscript 
denotes  summation  from  one  to  three,  and  Ft-  denotes  the  body  force  per  unit  volume  along 
the  xt  axis.  At  any  point  on  the  boundary  T  with  normal  v,  the  i  component  of  the  surface 
traction  per  unit  area  S\w\  can  be  written  from  the  condition  of  equilibrium  as 

S\v}  =  T«  V;  on  T, 

(4) 

where  Vy  is  the  cosine  of  the  angle  between  the  normal  v  and  the  x}  axis.  Substituting  (2)  into 
(3)  and  (4),  we  obtain 


(5) 

(6) 


S,yX0,y  +  2\ieyj  (SyA.0  j  +  2[LSjjJ)  +  F  ■  —  0, 

SM  =  Vy  [  5 ijke  +  2| xei}  -  (5(yX0"  +  2 \ie-j) } 

It  is  seen  that  -  (5,y^9,y  +  2(le,yy)  and  (5tJXQ  +  2|i0ey)  Vy  are  equivalent  to  F-t  and  S\v\  in  caus¬ 
ing  the  strain  field  and  are  here  denoted  by  F,  and  5,(v)  respectively,  giving 

5^0y  +  2peIyy  +  F1.  +  F;  =  O, 

5/(v)  +  s; (V)  =  Vy  (5yX0  +  2|iejy) 

Hence,  the  strain  distribution  in  a  body  with  inelastic  strain  under  external  load  is  the  same  as 
that  in  an  elastic  body  (no  inelastic  strain)  with  the  additional  equivalent  body  and  surface 
forces  Fi  and  S,(v).  This  reduces  the  solution  of  stress  field  of  a  body  with  known  inelastic 
strain  distribution  to  the  solution  of  an  identical  elastic  body  with  an  additional  set  of 
equivalent  body  and  surface  forces.  This  gives  the  same  results  as  Eshelby’s  famous  process 
of  imaginary  cutting,  relaxing,  restoring,  welding,  relaxing  in  his  noteworthy  paper  on  ellip¬ 
soidal  inclusions  [21], 


(7) 

(8) 
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If  there  is  thermal  strain  only  with  thermal  coefficient  of  expansion  a  and  temperature 
T,  we  can  write 

eV}  —  5y<x7;  0  =  =  3a T 

Then  the  equivalent  body  and  surface  forces  become 

F, 7  =  -(el  +  2\i)lTti;  =  (31  +  2p. 

This  is  the  well-known  Duhamel’s  analogy  between  temperature  gradient  and  the  body  force 
in  an  elastic  medium. 

Material  Defects  and  Initial  Stress  Field 

Imperfections  like  dislocations  exist  in  all  metals  and  cause  an  initial  stress  field.  Con¬ 
sider  a  most  favorably  oriented  crystal  at  a  free  surface.  For  a  slice  of  extrusion  to  form, 
positive  shear  deformation  must  occur  on  one  side  of  the  extrusion  and  negative  shear  on  the 
other.  The  initial  stress  field  x7  near  the  surface  (Fig.  1)  favorable  for  the  initiation  of  extru¬ 
sion  ”R'  is  one  with  positive  resolved  shear  stress  in  P  above  R  and  a  negative  one  in  Q 
below  R.  Such  an  initial  stress  field  can  be  provided  by  an  inelastic  strain  field  eaa  in  R. 
Consider  this  strain  to  be  linear  i.e.  is  constant.  This  causes  a  uniform  Fa  in  R  and 

hence  an  initial  positive  and  negative  resolved  shear  stress  x7  in  P  and  Q  respectively.  It  has 
been  pointed  out  by  Lin  and  Ito  [22]  that  this  eaa  can  be  produced  by  a  row  of  dislocation 
dipoles.  Recently  this  type  of  dislocation  dipoles  has  been  observed  under  microscope  in 
fatigue  specimens  as  indicated  by  Essmann  et  al.  [23]  and  Mughrabi  [24], 
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A  Gating  Mechanism 


Consider  a  most  favorably  oriented  crystal  at  a  free  surface  with  such  linear  e"aa  in  a 
polycrystal.  A  tensile  loading  causes  a  uniform  resolved  shear  stress  T4  in  the  crystal.  The 
signs  of  the  initial  shear  stress  and  the  applied  stress  are  of  the  same  sign  in  P  and  of  opposite 
sign  in  Q.  The  stress  x7  +  X4  in  P  is  largest  in  the  crystal.  Hence  P  reaches  the  critical  shear 
stress  xc  first  and  slides.  This  causes  a  residual  shear  stress  x*.  Due  to  the  continuity  of  the 
stress  field  t*,  slip  in  P  relieves  not  only  the  positive  shear  stress  in  P,  but  also  in  the  neigh¬ 
boring  region  including  Q.  This  keeps  the  positive  shear  stress  in  the  vicinity  from  reaching 
that  of  P  and  hence  from  sliding  in  the  tensile  loading.  This  relief  of  positive  shear  stress  in 
Q  is  the  same  as  the  increase  of  negative  shear  stress.  Hence  during  the  compressive  loading, 
Q  has  the  highest  negative  shear  stress  and  slides.  Again  this  negative  slip  causes  relief  of 
negative  shear  stress,  not  only  in  Q  but  also  in  P.  This  relief  of  negative  shear  stress  in  P  is 
the  same  as  increase  of  positive  shear  stress  causing  P  to  be  more  ready  to  slide  in  the  next 
tensile  loadings.  This  gives  the  natural  gating  mechanism  to  cause  alternate  sliding  in  P  and 
Q.  Hence  slips  in  P  and  Q  deepen  the  valley  and  raise  the  peak  monotonically  as  observed  in 
experiments.  These  microstress  fields  supply  a  natural  gating  mechanism  in  fatigue  bands. 
With  this  mechanism,  cross-slip  is  not  necessary  to  supply  of  irreversible  slip  as  indicated  by 
Mughrabi  [25], 

Consider  the  case  with  positive  e'aa  and  e^a.  Referring  to  Fig.  1,  the  equivalent 
force  Fa  is  pointing  outward  of  the  free  surface.  This  caused  an  extrusion.  As  the  extrusion 
grows,  the  thin  slice  R  increases  in  length.  This  elongation  cause  a  tensile  stress  xaa.  For  a 
f.c.c  crystal,  there  are  four  slip  planes.  This  xaa  causes  the  resolved  shear  stress  in  another 
slip  system.  To  simplify  the  numerical  calculation,  a  second  slip  system  is  assumed  to  have  a 
sliding  plane  with  a  normal  along  -axis  and  a  sliding  direction  along  x2-axis  or  vice-versa. 
The  resolved  shear  stress  x12  increases  with  x^.  After  x12  reaches  the  critical  shear  stress, 
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this  slip  system  slides.  The  slip  causes  a  plastic  strain  el2,  which  causes  a  plastic  tensile 
strain  eaa.  Now  the  total  inelastic  tensile  strain  is  equal  to  the  sum  e£,a  +  e'aa.  This  induces 
large  resolved  shear  stress  (initial  plus  residual)  of  opposite  signs  in  P  and  Q,  causes  a  higher 
rate  of  extrusion  growth,  and  explains  why  the  fatigue  initiation  growth  rate  is  higher  for  f.c.c 
metals  than  for  hexagonal  metals.  Some  calculations  of  this  secondary  slip  effect  were 
reported  by  Lin  and  Lin  [26] 

Mughrabi,  et.  al.  [25]  have  suggested  a  model  of  dislocation  dipoles  in  a  single  crystal 
as  shown  in  Fig.  2.  The  initial  inelastic  strain  due  to  the  interstitial  dipoles  is  pp  where 

m  denotes  the  number  of  these  dipoles.  Mughrabi  refers  to  mb  as  the  static  extrusion.  Cyclic 
loading  causes  extrusion  to  occur  on  both  the  right  and  left  ends  of  the  persistent  slip  band. 
Part  of  the  static  extrusion  is  lost  through  the  extrusion  process.  An  interesting  question 
which  was  raised  by  Mughrabi,  and  Essmann  et.  al.  [23]  was,  after  the  extrusion  has 
reached  the  amount  of  static  extrusion  ,  will  the  extrusion  growth  process  stop?.  Some 
details  of  the  answer  are  given  in  the  following  analysis. 

Now  we  consider  an  aluminum  polycrystal  of  fine  grain,  the  grain  size  is  small  as 
compared  to  the  total  volume  of  the  aggregate.  Hence  the  inelastic  strain  in  the  crystal  at  the 
free  surface  may  be  considered  to  be  occurring  in  a  semi-infinite  medium.  Referring  to  Fig. 
L  the  thickness  of  the  slices  P,  Q,  and  R  are  much  less  than  the  length  of  the  traces  of  the 
slip  lines  on  the  free  surfaces.  Consequently,  the  inelastic  strains  and  their  equivalent  forces 
are  assumed  to  be  constant  along  the  slip  line  direction,  and  the  semi-infinite  medium  is  taken 
to  be  under  plane  deformation. 

The  polycrystal  considered  is  pure  aluminum  subject  to  alternate  tension  and  compres¬ 
sion  %Q  along  x2-axis.  The  resolved  shear  stress  in  the  most  favorably  oriented  slip  system  is 
the  same  at  all  points  in  the  crystal. 


Fig.  2.  An  Interstitial  Dislocation  Dipole  Layer 
in  a  Single  Crystal 


10 


^af}  (i)  ~  ^2 


(9) 


The  plastic  strain  e ^  varies  from  point  to  point  in  the  slices  P  and  Q.  After  slip  occurs,  a  resi¬ 
dual  stress  field  is  produced.  From  the  plane  strain  solution  of  a  semi-infinite  medium,  the 
stress  field  is  written  as 


^ap  (*)  =  +  ^p  (*)  -2  G  £  C  (x  ,  ap;xn,Xp)  <?; 


aP- 


(10) 


where  xn  denotes  the  center  point  of  the  nth  grid  and  the  plastic  strain  of  the  nth  grid.  In 
the  active  grids  this  resolved  shear  stress  equals  the  critical  shear  stress  xc.  The  local  strain 
hardening  in  the  slip  bands  is  much  less  than  the  average  macroscopic  strain-hardening  of  the 
metal.  Neglecting  this  local  strain-hardening,  and  differentiating  (10)  with  respect  to  x0  yields 


£  xap(*) 

d  xQ 


V2  -2  G  £  C  (%  a(3;  xn  a P) 

n 


£  eap n 

d  x0 


(11) 


In  active  grids,  the  change  of  resolved  shear  stress  equals  the  change  of  critical  shear  stress. 


This  is  zero  due  to  the  neglect  of  local  strain  hardening.  Hence 
2  G  X  C(£  ,  ap;  xn,  ap)  5  eap-  = 

•  n  2  (12) 

There  are  as  many  unknown  5  <?apn  as  there  are  equations.  The  plastic  strain  increments 
8  <eap  n  in  sliding  grids  for  an  incremental  applied  stress  5  x0  can  be  determined.  Thin  slices 
®  P,Q,  R  as  shown  in  Fig.  1  are  taken  to  be  0.05  um  each  along  x2-axis.  The  center  distance 

between  P,Q  is  0.10  um  and  the  linear  dimension  of  the  crystal  is  take  to  be  50  um.  For 
numerical  calculation,  the  crystal  is  divided  into  1010  parallelogram  grids  oriented  45°  to  the 

•  free  surface  with  10  equal  spacings  along  jcraxis  and  101  rows  along  *2-axis.  The  plastic 
strain  in  each  grid  is  assumed  to  be  uniform.  The  stress  field  caused  by  a  uniform  plastic 
strain  in  a  grid  centered  at  (x)  was  calculated.  The  relief  of  the  resolved  shear  stress  at 
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C x )  due  to  this  eap  is  called  residual  shear  stress  expressed  as 


*£p  ( *.)  =  2GC  (x  ,ap;  x„,aP) 


(13) 


where  2GC  (jr,  aP;  x„,  ap)  is  the  resolved  shear  stress  xap  at  x  due  to  unit  plastic  strains  ea$ 
at  xn.  This  is  called  the  resolved  shear  stress  influence  coefficient.  The  detail  of  the  calcula¬ 
tion  of  this  influence  coefficient  is  shown  in  by  Lin,  [27],  and  Lin  &  Lin  [26]. 


Effect  of  Primary  Slip  on  Stress  in  Secondary  Slip 


One  set  of  uniform  initial  shear  stress  of  24.4  KN/m 2  (3.53  p.s.i)  in  P  and  -24.4  KN/m2 
(-3.52  p.s.i.)  in  Q,  and  applied  tensile  and  compressive  stresses  of  +690.0  KN/m 2  (+  100  p.s.i.)  * 

was  calculated  by  Lin  and  Lin,  [28].  After  200  cycles,  the  plastic  strain  distributed  in  P,Q 
were  obtained  and  are  shown  in  Fig.  3.  This  elastic  strain  is  much  less  than  these  plastic 
ones.  Hence  these  plastic  strains  represent  the  total  shear  strain  in  P,Q.  Let  the  thickness  of 
P,Q,  R  be  t,  then  t  =  0.5/V2  \xm.  The  displacement  in  R  along  a-direction 

2(eap)p,g  r  =  2(ea p)p,2  r 

04)  • 

This  is  similar  to  shearlag  problems  in  thin  wall  airplanes  structures. 


The  residual  stress  (Tpp  (x)  may  be  expressed  as 
CTpp(i)  =  2GC  (x  ,PP;  x*.  ap)e"Pi> 

Using  these  influence  coefficients,  C  (x,  PP;  x  „,  ap)  it  was  found  Opp  (x)  in  R  is  very  small. 
*^PP  ©  (*) 
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Fig.  3.  Plastic  Shear  Strain  Distributions 
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The  resolved  shear  stress  x12  in  the  secondary  slip  plane  in  R  is  then 


xl2  - 


'aa 


(20) 


Effect  of  Secondary  Slip  on  Primary  Slip 

When  x12  reaches  tF  in  some  part  of  R,  e[2  occurs  in  that  region.  This  e[2  has  a  strain 

//  f 

component  eaa  which  adds  to  the  in  causing  the  difference  in  resolved  shear  stress  xap  in 
P  and  Q.  This  increases  the  rate  of  sliding  in  P,Q  thus  accelerating  the  extrusion  growth  at 
the  free  surface  as  reported  by  Lin  and  Lin  [26].  This  additional  eaa  caused  by  el2  in  the 
second  slip  system  in  R  causes  the  extra  length  of  slice  R  to  be  more  than  the  initial  static 
extrusion.  This  is  one  of  the  main  reasons  why  the  extrusion  can  proceed  on  beyond  the 
static  extrusion  as  observed  by  Essmann  et  al.  [23]. 

The  distribution  of  e12  in  R  depends  on  the  variation  of  slope  in  e'a^  in  P  and  Q.  The 
latter  depends  on  the  distribution  of  the  initial  resolved  shear  stress  in  P,Q.  To  simplify  the 
numerical  calculation,  it  was  assumed  that  only  positive  slip  c”p  occurred  in  P,  only  negative 
slip  -eap  occurred  in  Q  and  only  e12  occurred  in  R.  Three  cases  of  initial  resolved  shear 
stress  were  calculated.  They  are  as  follows. 


1)  The  initial  resolved  shear  stress  varies  linearly  from  the  interior  grain  boundary  to  the  free 
surface 

t7  =  +34.5  -^  ( 1-4)  »+  in  p »  and  -  in  Q. 

m 2  d 

The  plastic  strains  in  P,Q  are  about  the  same.  This  distributions  of  primary  slip  in  P,Q  with 
and  without  the  consideration  of  the  secondary  slip  in  R  are  shown  in  Fig.  4.  The  distribu- 
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Plastic  strain 


Fig.  4.  Plastic  Shear  Strain  e%  in  P,Q 
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tion  of  the  secondary  slip  in  R  is  shown  in  Fig.  5.  It  is  to  be  noted  that  e[2  curve  in  R  is  very 
close  to  the  slope  of  ea p  in  P,Q.  The  equivalent  body  force  caused  by  e'[2  is  proportional  to 
the  slope  of  e[2  curve  in  R  and  hence  to  be  the  curvature  of  in  P,  Q. 

2)  The  initial  resolved  shear  stress  in  P,Q  varies  linearly  from  a  maximum  value  of  3.45 
KN/m2  (5  p.s.i)  at  the  free  surface  to  13.5  KN/m 2  (2  p.s.i)  at  the  grain  boundary.  The  distribu¬ 
tions  of  the  primary  slip  in  P,Q  with  and  without  the  consideration  of  the  secondary  slip  in  R 
are  shown  in  Fig.  6.  The  variation  of  the  secondary  slip  in  R  is  shown  in  Fig.  7. 

3)  The  initial  resolved  shear  stress  in  P,Q  varies  parabolically  from  zero  (with  zero  slope)  at 
the  grain  boundary  to  a  maximum  value  of  34.5  KN/m 2  (5.  p.s.i.)  at  the  free  surface.  The  dis¬ 
tribution  of  the  primary  slip  in  P,Q  with  and  without  the  consideration  of  the  secondary  slip  in 
R  and  shown  in  Fig.  8.  The  variation  of  secondary  slip  in  R  is  shown  in  Fig.  9. 

The  growth  of  plastic  strain  e'a^  at  the  free  surface  with  cycles  of  loading  for  the  three 
cases  of  initial  stresses,  with  and  without  secondary  slip,  is  shown  in  Fig.  10. 

Extrapolation  Proceedure 

The  calculation  of  the  secondary  slip  may  require  computations  for  a  large  number  of 
cycles  of  loadings.  To  reduce  the  computation,  an  extrapolation  proceedure  was  developed. 
At  the  end  of  the  compression  loading  of  the  Nth  cycle,  the  incremental  plastic  strain  in  P,  Q 
and  R  in  the  Nth  cycle  were  calculated  and  denoted  by  To  extrapolate  to  the  (N+n)th 
cycle,  the  plastic  strains  in  P  were  assumed  to  increase  by  n  x  Since  at  the  end  of  the 
(N+n)th  cycle,  active  grids  in  Q  and  R  just  finish  sliding  and  hence  have  stresses  equal  to  the 
critical  shear  stress  xc.  We  have 

T'4  +  x7  +  T*  =  xc  in  Q  and  R 

(17) 

This  condition  was  used  to  find  the  incremental  plastic  strains  in  Q  and  R.  The  number  of 
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Fig.  5.  Plastic  Shear  Strain  ej^in  R 
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Fig.  6.  Plastic  Strain  Distribution  e^p  in  P.Q 
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Fig.  7.  Plastic  Strain  Distribution  e^in  R 
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Fig.  8  -  Plastic  Distribution  e12  in  R. 
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Fig.  9  -  Plastic  Strain  Distribution  ea^  in  P,Q. 
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Fig.  10  -  Growth  of  Plastic  Strain  at  Free  Surface 
With  and  Without  Secondary  Slip 


equations  in  the  form  (17)  equals  the  number  of  unknowns  of  (Se)’s  in  Q  and  R.  These 

_  rr  'i 

be  ’s  were  thus  determined.  The  number  of  cycles,  n,  for  each  extrapolation  was  10  cycles. 
The  results  for  the  first  extrapolation  were  compared  with  those  calculated  cycle  by  cycle  and 
the  agreement  was  found  to  be  good. 

Some  Related  Metallurigical  Observations 

A  number  of  metallurgical  observations  supporting  the  micromechanic  theory  of 
fatigue  crack  initiation  have  been  shown  by  Lin  (1977,  1981,  1984).  Some  metallurgical 
observations  by  Mughrabi  (1983),  related  to  the  present  study,  are  given  as  follows: 

1.  The  persistent  slip  bands  PSB  in  pure  f.c.c.  metals  consist  of  more  or  less  equally 
spaced  dislocation  walls,  which  are  densely  packed  with  edge  dislocation  dipoles  and 
which  are  separated  by  so  called  channels  of  low  dislocation  density.  These  disloca¬ 
tion  dipoles  can  cause  different  initial  shear  stresses  of  opposite  signs  in  P,Q.  Three 
sets  of  initial  stresses  are  here  shown. 

2.  Single  copper  crystal  fatigue  tests  show  that  the  extrusion  height  has  exceeded  the 
static  extrusion  by  a  large  amount.  As  indicated  by  Mughrabi  these  observations  do 
not  support  the  simple  theory  of  the  termination  of  extrusion  growth  after  extrusion 
height  reaches  the  static  extrusion.  The  rate  of  growth  does  not  depend  markedly  on 
temperature,  and  the  absence  of  this  dependence  suggests  that  some  micromechanics 
process,  as  proposed  here,  seems  to  be  the  cause  of  extrusion  growth  beyond  the  static 
extrusion. 

3.  TEM  (Tramission  Electron  Microscope)  observations  show  that  late  in  the  fatigue  test, 
increasing  secondary  slip  activity  in  the  PSB’s  gradually  transforms  the  primary  edge 
dislocation  PSB  walls  into  a  cell  structure.  This  shows  that  secondary  slip  occurs. 


4. 


The  extrusion  height  in  a  copper  polycrystal  of  a  grain  size  of  25  \1  was  estimated  to 
reach  O.lp  to  0.2|i,  indicated  by  Mughrabi  (1983).  For  a  0.035)1  thickness  of  slices 
P,Q,  used  in  the  present  calculation,  a  of  about  2.8  will  yield  an  extrusion  of 
height  of  O.lp.  The  calculated  eap  at  the  free  surface  about  1.5.  It  is  expected  that 
increase  of  loading  cycles  can  yield  a  of  2.8. 

Conclusions 

Three  thin  adjacent  slices  P,R,Q,  in  a  crystal  at  the  free  surface  with  initial  resolved 
shear  stress  of  opposite  signs  in  P,Q  were  assumed  in  a  polycrystal.  This  polycrystal  was 
subjected  to  cyclic  tension  and  compression  loading.  Primary  slip  in  P,Q  induces  a  tensile 
stress  and  a  resolves  shear  stress  in  a  secondary  slip  system  in  R.  When  the  shear  stress 
reaches  critical,  this  secondary  slip  system  slides  and  gives  plastic  shear  strain  e[2  which  has  a 
strain  component  eaa.  The  eaa  has  the  same  effect  in  causing  the  initial  stress  of  opposite 
signs  in  P,Q  as  the  initial  inelastic  tensile  strain  which  may  be  caused  by  interstitial  dislo¬ 
cation  dipoles,  shown  by  Lin  Ito  [22]  and  Lin  and  Lin  [18].  Hence,  the  secondary  slip  causes 
the  primary  slip  to  proceed  further.  eaa  contributes  to  the  extrusion  growth  beyond  the  intitial 
static  extrusion.  The  same  reasoning  can  be  applied  to  intrusion  growth  by  an  array  of 
vacancy  dipoles  instead  of  the  interstitial  ones.  After  the  intrusion  is  developed  to  a  certain 
size,  the  stress  concentration  will  contribute  to  the  crack  initiation. 
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